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We investigate the superconducting mechanism and the transition temperature of heavy 
fermion superconductor UPd2Al3 on the basis of a single band two-dimensional Hubbard 
model on triangular lattice, which represents the most heavy band of UPd2Al3. Both normal 
and anomalous self-energies are calculated up to third order with respect to the Coulomb 
repulsion U between itinerant electrons. The superconducting transition temperature is ob- 
tained by solving the Eliashberg's equation. Reasonable transition temperatures are obtained 
for moderately large U. It is found that the momentum and frequency dependence of spin 
fluctuations given by RPA-like terms gives rise to the d-wave pairing state, while the ver- 
tex correction terms are important for obtaining reasonable transition temperatures. These 
results seem to show that the superconductivity in UPd2Al3 can be explained by the per- 
turbation theory with respect to U. 

KEYWORDS: heavy fermion systems, electron correlation, UPd2Al3, superconducting transi- 
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1. Introduction 

1.1 Introduction to UPd2Al3 

UPd2Al3 crystallizes in the hexagonal PrNi2Al3-type structure and it was found to be a 
superconductor with Tc =2K. The superconducting state coexists with antiferromagnetic order 
below the Neel temperature Tn=14.5K.^ The magnetic moments are aligned ferromagnetically 
in the hexagonal basal plane and coupled antiferromagnetically between planes with wave 
vector k = [0,0,0.5].^ 

For T < Tn, the electrons in UPd2Al3 behave as if they were separated into two rather 
independent subsystems of 5f character.^' ^ One subsystem is a localized system which is 
responsible for antiferromagnetic order. We call this system "localized subsystem" . The second 
subsystem is the itinerant one characterized as a heavy electron system and responsible for the 
superconductivity. We call the latter "itinerant subsystem" . For Tc < T <C Tn in the normal 
state, the electronic specific- heat C can be approximated by C ~ 7T where 7 = 145m J/K-^.^ 
For Tc < T <4K in the normal state, the electrical resistivity curve follows the dependence 
P = Po + where po = 2.3/xr2cm and A ~ 10~^/ir2cm/K^.^ The behavior observed in the 
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electronic specific-hcat and the electrical resistivity indicates that the itinerant electrons in 
UPd2Al3 compose a Fermi liquid. 

The calculated Fermi surface in the antiferromagnetic state has four principal sheets 
of predominately /—electron character. These are corroborated by de Haas-van Alphen 
measurements.® All bands arc doubly degenerated, because the antiferromagnetic structure 
is conserved under a time-reversal and inversion-symmetry operator with the non-primitive 
translation vector v = (0,0,0.5). The Fermi sheets are not of hexagonal symmetry, reflecting 
the antiferromagnetic structure. Roughly speaking, two Fermi sheets have two-dimensional 
topology, and the remaining two Fermi sheets have the dispersion along the crystallographic 
c— axis, although there exist some discrepancies between the results of ref.7 and ref.8 . In ref.7, 
two Fermi sheets which have two-dimensional topology are called , "party hat" and "cylinder", 
and the remaining two are called "cigar" and "eggs". According to the band calculation and 
de Haas-van Alphen measurements,® the band which has the largest effective mass is "party 
hat" . 

As for the superconducting property, the specific heat and the nuclear spin-lattice relax- 
ation rate Ti show the T^- dependence without the coherence peak just below Tc-^^^^"^ The 
Knight shifts decrease below Tc-^^ The thermal conductivity k dose not show any anomaly at 
Tc and shows the behavior below IK.^ These facts suggest that the superconductivity is in 
principle of a d-wave type characterized by vanishing gaps on the lines on the Fermi surface, 
although there exist some discrepancies between the results of the specific heat and NMR. 
Inelastic neutron scattering i-^^i*^ arid tunneling experiments-'^^' below the superconducting 
transition temperature have been performed on this material. The results of these experiments 
indicate that superconductivity and magnetic excitations of antifcrromagnetically ordered mo- 
ments in UPd2Al3 couple each other. 

1.2 Motivation 

Itinerant electrons in heavy fermion systems compose a Fermi liquid whose effective mass 
is enhanced by 10^ ~ 10^ times as large as that of free electron. The electron correlation 
is essential in giving the large electron mass in heavy electron systems. From this point of 
view, we can expect that superconductivities in heavy fermion systems are derived also from 
the electron correlation through the momentum and frequency dependence of the effective 
interaction between electrons. In fact, in other strongly correlated electron systems such as 
high-Tc cuprates,^^ organic superconductors^^ and spin-triplet superconductor Sr2Ru04,^^ it 
has been recognized that the electron correlation plays an important role in realizing their 
superconductivities, by treating the electron correlation by perturbation theory based on 
Fermi liquid theory. 

As described above, recent neutron scattering experiments done below Tc revealed a strong 
relationship between the superconductivity and magnetic excitations of antiferromagnetically 
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ordered moments in UPd2Al3. Indicated by these experimental faets, several groups^'^"^^ have 
proposed a new mechanism mediated by spin wave in the ordered antiferromagnetic state. 
This mechanism is called "magnetic exciton" mechanism. Magnetic excitons are a kind of spin 
waves that ripple through a magnetically ordered medium, but they require a certain amount 
of energy to be excited. They argued that these magnetic excitons may produce an effective 
interaction between the itinerant electrons, and so be responsible for superconductivity in 
a manner analogous to that played by phonons in conventional superconductors. Therefore, 
for "magnetic exciton" mechanism, the coupling between the "localized subsystem" and the 
"itinerant subsystem" is essential to deriving the superconductivity. 

The proposal seems to be interesting as a new mechanism. However, it is not clear whether 
this mechanism can actually realize the superconductivity or not, since there is no quantitative 
calculation of the critical temperature Tc. To the present authors, it seems to be natural that 
the largest Fermi sheet with large electron mass itself leads to the superconductivity by a 
mechanism due to the effective interaction between electrons in the "itinerant subsystem". 

In this paper, we investigate the possibility that the superconductivity in UPd2Al3 is 
derived from the electron correlation effects existing in "itinerant subsystem" , by the pertur- 
bation theory with respect to Coulomb repulsion U. We show that the superconductivity in 
UPd2Al3 can be explained by the usual spin fluctuation mechanism in the "itinerant subsys- 
tem" . This explanation will be given on the basis of detailed calculation where the dependence 
on the parameters inherent in UPd2Al3 are carefully studied. We concentrate our attention 
on the calculation of the superconducting transition temperature Tc, so we do not investigate 
the relationship between superconductivity and magnetic excitation of antiferromagnetically 
ordered moments in UPd2Al3 far below Tc, in this paper. 

The organization of this paper is as follows. In §2 formulation is presented, including 
model Hamiltonian, perturbation expansion terms and Eliashberg's equation. In §3 calculation 
results of Tc and other various quantities are shown. Finally in §4 discussions and conclusion 
are presented. 

2. Formulation 

2. 1 Model and Hamiltonian 

Taking the importance of the correlation between electrons in the "itinerant subsystem" 
into account in realizing the d-wave superconductivity, we consider only a single band of 
UPd2Al3 possessing the largest effective mass. According to the band calculation,^ the band 
of UPd2Al3 possessing the largest effective mass (which is the so called "party hat") has 
the property of quasi-two-dimensional. Therefore we represent the band as a simple effective 
two-dimensional Hubbard model on the triangular lattice. We set x-y plane in the hexagonal 
basal plane, and set y-axis parallel to antiferromagnetic moments. We consider only the nearest 
neighbor hopping integrals and we assume for the ^/-direction hopping integral a different value 
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tm from other hopping integrals t, because the superconductivity of UPd2Al3 is reahzed in the 
antiferromagnetic statc^(see Fig. 1). So we consider that the effect of the antiferromagnetic 
order(the effect of the "locahzed subsystem") is included in the difference between and t, 
and we determine the values so as to reproduce the considered Fermi sheet which is obtained 
by the band calculation and is not of hexagonal symmetry, reflecting the antiferromagnetic 
structure ( "localized subsystem" ) . 

We rescale length, energy, temperature, time by a,t,-^,j respectively (where a,kB,h 
are the lattice constant of the hexagonal basal plane, Boltzmann constant, Planck constant 
divided by 2-k respectively), we write our model Hamiltonian as follows: 

H = Ho + Hi, (1) 

Ho = J2{e{k)-,x)al^ak,, (2) 

k;,cr 

a/S 1 

e(k) = -4:Cos{—kx) cos{-ky) - 2tmC0s{ky), (3) 
^ ^ XI X '^ki+k2,k3+k4aki^a^^,ak30-'ak4(T, (4) 

tr^ir' kik2k3k4 

where a|^^(ako-) is the creation (annihilation) operator for the electron with momentum k 
and spin index a; tm and fi arc the y-direction hopping integral and the chemical potential, 
respectively. The sum over k indicates taking summation over a primitive cell of the inverse 
lattice. We set this primitive cell of the inverse lattice as shown in Fig. ?? for the convenience 
of our numerical calculation. 



2.1.1 model parameters 

Our model parameters are temperature T, j/-direction hopping integral t^, the Coulomb 
repulsion U and the electron number n per one spin site. The bare bandwidth Ae(tin) of our 
model is 

. ,^ . /8 (0<t^<0.5). 
[4 + 4t^ + l/U (tm>0.5). 
According to the band calculation and the de Haas-van Alphen experiment,^the parameter 
region of our model which reproduces the Fermi sheet is given by 0.50 < < 0.75 and 
0.51 <n< 0.572. In this parameter region, the range of the bare bandwidth is 8 < Ae(tin) < 
8.333 • • • . We consider that t-ai = 0.75, n = 0.572 are the best fitting values of parameters which 
well reproduce the Fermi sheet we consider. In Fig. 3, the bare Fermi surface calculated by our 
model (the model parameters are tm = 0.75, n = 0.572, U = and T = 0) is presented. From 
this figure we can see that the considered Fermi sheet determined by the band calculation^ is 
well reproduced within our effective model. 
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The superconducting transition temperature determined by experiments is about 
2K.^Then the rescaled superconducting transition temperature Tc is Tc — 1.4 x 10"'^, when 
we set leV as the value of bandwidth obtained by the band calculation.^'^ 

2.2 Green's functions 

2.2.1 Bare Green's functions 

The bare Green's function is the following, 

= : \- , (6) 

where e„ = (2n + l)TrT{n : integer) is the fermion-Matsubara frequency and short notation 
k = (k, e„) is adopted. We consider the Hartree term is included in the chemical potential. 

2. 2. 2 Dressed normal Green's functions 

Next we consider the dressed normal Green's function. When we consider the situation near 
the superconducting transition temperature, the Dyson-Gorkov's equation can be linearized. 
Therefore the dressed normal Green's function is obtained from the bare Green's function 
with only the normal self-energy correction. We expand the normal self-energy up to third 
order with respect to U ;the diagrams are shown in Fig. 4. 

Then we obtain the normal self-energy as follows, 

^-(k) = ^ J2P^Xo{k - k') + U\l{k - k') + U^4>lik + k')]Go{k'), (7) 

k' 

where Xo(' " ' ) <;i!>o(- • • ) are given respectively as 

X0(q,'^m) = -J^^GQ{]s.,en)GQ{(\ + ]^,Um + en), (8) 
k,n 

^o(q,^^m) = -j;^^GQ{)^,en)Go{(i-\<.,LOra - en)- (9) 
k,n 

Here oJm = 2rmrT(m : integer) is the boson-Matsubara frequency. The quantity Xoi^^'^m) 
has the physical meaning of the bare susceptibility and expresses spin fluctuations in the 
system. More over, the bare susceptibility plays an important role in the calculation of Tc, 
that is, it determines the magnitude and the spatial and temporal variation of the effective 
interaction between electrons in the "itinerant subsystem" , through the higher order terms in 
U. (See the first term of right hand side of the equation 2.13.) 

Note that the Hartree term has been already included in the chemical potential and the 
constant terms which have not been included in the Hartree term are included in the chemical 
potential shift when we fix the particle number. 

Then the dressed normal Green's function is 

^^^^ ^ i6„-(e(k)-/x-<5/. + S„(fc))' ^^^^ 
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where S/j, is determined so that the following equation is satisfied. 

^=^T.Gik) = ^Y.Go{k). (11) 

k k 

We expand the above equation up to the third order of the interaction with regard to 611 — 
T,n{k), we obtained 6n as 

^~ Xo(0,0) • ^^'^ 

2.2.3 Anomalous self-energy and effective interaction 

When wc consider the situation near the superconducting transition tcmpcratiuc, the 
anomalous sclf-cncrgy Tia{k) is represented by the anomalous Green's function F{k) and 
the (normal)cffective interaction. We expand the effective interaction up to third order with 
respect to U as shown in Fig. 5. 

Then we obtain the anomalous self-energy as follows; 

^«(^) = -'^Y.\!^ + U\o{k + k') + 2U\l{k + k')\F{k') 

k' 

- E Go{k')[xo{k + k') - Mk + k')]Go{k + k' - k")F{k") 

k'k" 

- E Go{k')[xo{-k + A;') - 4>^{-k + fc')]Go(-A; + k' - k")F{k"). (13) 

k'k" 

2.3 Eliashberg's equation 

Prom the linearized Dyson-Gorkov equation, we obtain the anomalous Green's function 
as follows; 

F{k) = \G{k)fj:a{k). (14) 
Then the Eliashberg's equation is given by 

^«(^) = -^J2iU + U\o{k + k') + 2U^xlik + k')]\Gik')\^^a{k') 

k' 

- E Goik')[xo{k + k') - Mk + k')]Go{k + k' - k")\G{k")\^^a{k") 

k'k" 

- E Go{k')[xo{-k + k') - M-k + k')]Go{-k + k' - k")\G{k")\''^a{k"). 

k'k" 

(15) 

We consider that the system is superconducting state when the eigen value of this equation 
is 1. 
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3. Calculation Results 

3.1 Details of the numerical calculation 

To solve the Eliahberg's equation by using the power method algorithm, we have to cal- 
culate the summation over the momentum and the frequency space. Since all summations are 
in the convolution forms, we can carry out them by using the algorithm of the Fast Fourier 
Transformation. For the frequency, irrespective of the temperature, we have 1024 Matsubara 
frequencies. Therefore we calculate throughout in the temperature region T > Tn^ , where 
Tiijn is the lower limit temperature for reliable numerical calculation, which is estimated about 
2.0xl0~^(> Ae(tm)/(27rxl024) ~ 1.3 x 10"^); we divide a primitive cell into 128x128 meshes. 

We have carried out analytically continuing procedure by using Fade approximation. 

3. 2 Dependence of Tc on U, tm and n 

To solve the Eliashberg's equation, we set the initial gap function((ia;y-symmetry) as fol- 
lows. 

\/3 1 

T,a{k) oc sin -^k^ sin -ky. (16) 

Notice that the "itinerant subsystem" has no hexagonal symmetry due to the antiferromag- 
netic structure( "localized subsystem"). The calculated gap functions show the node at k^ = 
and ky = and changes the sign across the node for all parameters. The symmetry of Cooper 
pair is dxy 

We calculate Tc around the best fitting values of parameters (tm = 0.75, n = 0.572) because 
the best fitting values of parameters have some arbitrariness. 

The dependence of Tc on U, n and tm are shown in Figs. ??~ 9. Prom these results, we can 
point out the following facts. For large U higher Tc are obtained commonly for all parameters. 
When we fix the electron number n per one spin site and the Coulomb repulsion U, for large 
fm, the system get close to real triangular lattice {t^ = 1) for electrons in the "itinerant 
subsystem" and at the same time Fermi level goes away from the van Hove singularity (see 
Fig. 14 shown later), then lower Tc are obtained. Contrary, when we fix the y— direction 
hopping integral and the Coulomb repulsion U, for large n, the system becomes away from 
the half-filling state (n = 0.5) and at the same time Fermi level gets close to the van Hove 
singularity (see Fig. 15 shown later), then higher Tc are obtained. 

3.3 Vertex corrections 

To examine how the vertex corrections influence Tc, we calculate Tc by including only 
RPA-like diagrams of anomalous self-energies up to third order, in other words, without the 
vertex corrections and compare obtained Tc with Tc calculated by including full diagrams 
of anomalous self-energies up to third order(Fig. 10). From this figure, we can see that Tc 
calculated by including only RPA-like diagrams of anomalous self-energies up to third order 
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is higher than calculated by including full diagrams of them. This figure shows that the 
main origin of the d-wave superconductivity is the momentum and frequency dependence of 
spin fluctuations given by the RPA-like terms; spin fluctuations exist in our nested Fermi 
sheet, and the Coulomb interaction among electrons in "itinerant subsystem". The vertex 
corrections reduce by one order of magnitude. So the vertex corrections is important for 
obtaining reasonable T^. 

3.4 Behavior of xo{q) 

The calculated results of the static bare susceptibility are shown in Fig. 11~ 13 for various 
value of tm and n. Prom these figures, we point out the following facts. For all parameters, 
Xo(q;0) have sufficiently strong momentum dependence, although there exist no prominent 
peak. When we fix the electron number n per one spin site and increase the parameter tm, 
the system approaches the real triangular lattice for electrons in the "itinerant subsystem". 
In this case the peak and the momentum dependence of Xo(q)O) is slightly suppressed since 
the antiferromagnetic spin fluctuations are suppressed by the frustration effect in triangular 
lattice. Contrary, when we fix the y— direction hopping integral tm and increase the electron 
number n per one spin site, the system become away from the half-filling state and the peak 
and the momentum dependence of XoCQjO) is slightly enhanced. The reason is that Fermi 
surface possesses some nesting properties with increasing n from half-filling case, where the 
Fermi surface is similar to a circle without nesting properties. 

3. 5 Density of states 

The density of states(DOS) is given by 

Pi^) = -^T.^^G''(K^), (17) 

k 



where 



G^{k,uj) = G(k,e, 



We show the tm and n- dependence of DOS in Figs. 14 and 15. From insets in these figures, 
we can see that the position of the van Hove singularity shifts upward from the Fermi level 
with increasing tm and decreasing n. This departure of the van Hove singularity from the 
Fermi level reduces the superconducting transition temperature Tc- 

3.6 Self- energy 

The self-energy is given by 

The real part and the imaginary part of the self-energy at Fermi momentum are shown in 
Fig. 16 and Fig. 17 respectively. The w-dependence of both parts near a; = are respectively 
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given by RcS^(kj,a;) oc —lo and ImS^(kj,Lj) oc — o;^. This behavior is the same as that for 
the usual Fermi hquid. As U increases, the slope of ReS^(kj,ci;) at cj = beeomes steeper 
and the coefficient of the oj^-term in ImS^(kj, cj) become larger. This indicates that the mass 
and the damping rate of the quasi-particle become larger as U increases. These results are the 
typical Fermi liquid ones. 

4. Summary, Discussion and Conclusion 

In this paper, taking the importance of the correlation effect between electrons in "itin- 
erant subsystem" into account in realizing the superconductivity, we have considered a single 
band of UPd2Al3 whose effective mass is the largest. We have represented the band with an 
effective two-dimensional Hubbard model on triangular lattice and calculated on the basis 
of third order perturbation theory. Reasonable transition temperatures have been obtained 
for moderately large U. We also have calculated Tc by including only RPA-like diagrams of 
anomalous self-energies up to third order, in other words, without the vertex corrections and 
compared it with Tc calculated by including full diagrams of anomalous self-energies up to 
third order. We point out that the main origin of the superconductivity can be ascribed to the 
momentum and frequency dependence of the spin fluctuations and the correlation between 
electrons in the "itinerant subsystem", although the vertex correction terms are important 
for reducing and obtaining reasonable transition temperatures. We have calculated DOS 
and the normal self-energy. The obtained behavior of both quantities is the expected one for 
the Fermi liquid. 

Now, we briefly discuss the theories suggesting that the superconductivity in UPd2Al3 
is derived from the mechanism of exchanging "magnetic excitons". Huth et al.^ predicted 
that line nodes are at the rim of the magnetic Brillouin zone and on the "cigar" Fermi sheet. 
Their theory is based on the assumption of the pairing interaction mediated by the exchang- 
ing "magnetic excitons" and the band structure calculated by K.Knopfie et alJ Miyake and 
Sato^^ also predicted that the line node is on the plane very close to the zone boundary of the 
folded Brillouin zone in the antiferromagnetically ordered state on the basis of the itinerant- 
localized duality model; their theory is also based on the observed behavior of the dynamical 
magnetic susceptibility and the assumption of the pairing interaction induced by the exchang- 
ing "magnetic excitons" but without using any information about the Fermi sheets. According 
to these analysis,^^'^^ Fermi sheets which have the dispersion along the crystallographic c-axis 
and cross the zone boundary plane of the folded Brillouin zone in the antiferromagnetically 
ordered state are necessary to possess the line nodes actually. On the other hand, according 
to the band calculation performed by Inada et al.,^ the "cigar" Fermi sheet dose not cross 
the zone boundary plane of the folded Brillouin zone in the antiferromagnetically ordered 
state, while the "cigar" Fermi sheet calculated by Knopfie et alJ crosses this plane. Moreover, 
Inada et aZ.^'^^ have tried to determine the position of the line node on the actual Fermi 
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surface by the de Haas-van Alphen experiment. Their conclusion is that any line node in 
the anisotropic energy gap does not exist on the "cigar" Fermi sheet and/or it is difficult to 
distinguish experimentally the line node from the gapless state created by the pair-breaking 
effects by the applied magnetic field. Therefore, we now consider that the basis of "magnetic 
exciton" mechanism has not been confirmed. 

These theories described above are indicated by the existence of the coupling between the 
superconductivity and the magnetic excitations of antiferromagnetically ordered moments in 
UPd2Al3 below T^. Generally speaking, the coupling between superconductivity and mag- 
netism exists more or less. Therefore it doesn't always follow that superconductivity is mainly 
derived from the magnetic mechanism by the reason why there exists the coupling between 
superconductivity and magnetism. This is understood from the fact that the coupling between 
superconductivity and spin fluctuations or phonons is detected in high-Tc superconductors but 
the superconductivities in high-Tc are not always mainly derived from these spin fluctuations 
or phonons. 

On the other hand, above Tc, the behavior of the electronic specific- heat and the electrical 
resistivity are the typical Fermi liquid ones and anomalous behavior of these quantities is not 
detected in UPd2Al3 ,for example, such a,s p — ccT vn high Tc superconductor. Moreover the 
magnetic transition temperature (Tn = 14K) is sufficiently higher than the superconducting 
transition temperature (Tc = 2K). Therefore we consider that the coupling between local 
moments ( "localized subsystem" ) and electrons in the "itinerant subsystem" is not so strong, 
while Sato et al. assume that this coupling is very strong in ref.26. 

Therefore we considered that the electron correlation in the "itinerant subsystem" also 
plays an important role to derive the superconductivity in UPd2Al3. Quasi- two-dimensional 
Fermi sheet so called "party hat" is most favorable to superconductivity, because of the 
most strong electron correlation in the "itinerant subsystem". Our proposition is that this 
sheet possessing large area and electron mass is essential to derive the superconductivity in 
UPd2Al3, while the Fermi sheets called "cigar" and "eggs" playing an essential role in realizing 
the superconductivity mediated by "magnetic exciton" seem to be too small to induce the 
superconductivity in the total system. 

In the last, we also point out that results of tunneling experiments and our proposition 
don't contradict. Tunneling spectra have been also measured in this material by the same 
group. ^^'^'^ Their conclusion is that an energy gap along the c-axis in UPd2Al3 is well resolved. 
In the multi band system such an UPd2Al3, the superconducting transition, in principle, 
simultaneously occurs in all Fermi sheets. Therefore the energy gaps open in all Fermi sheets 
below Tc, then an energy gap along the c-axis exists in UPd2Al3 since UPd2Al3 is not a perfect 
two-dimensional system. 

In conclusion we present a new mechanism that the superconductivity in UPd2Al3 can 
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be also derived from the correlation between electrons in the "itinerant subsystem" and the 
symmetry of the pairing state is dxy ■ In our model, the main origin of the superconductivity is 
the momentum and frequency dependence of the spin fluctuations and the Coulomb interaction 
between electrons in the "itinerant subsystem" . This momentum and frequency dependence of 
the spin fluctuations stems from the shape of our Fermi sheet which undergoes the symmetry- 
breakdown due to the antiferromagnetic order ("localized system") and then possesses nesting 
properties. 
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Fig. 1. Schematic figure of the hexagonal basal plane of UPd2Al3. Circles in this figure represent U 
atoms. 




kx 



Fig. 2. Our primitive cell of the inverse lattice. The regions enclosed by the solid and dashed lines 
are our primitive cell and the first Brillouin zone, respectively. 



13/?? 



J. Phys. Soc. Jpn. 



Full Paper 




Fig. 3. The bare Fermi surface calculated by our model. The region enclosed by the dashed line is 
the first Brillouin zone. See also Fig. 5(b) in ref. 8. 





Fig. 4. The Feynman diagrams of the normal self-energy up to third order. The solid and dashed 
lines correspond to the bare Green's function and the interaction, respectively. 
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Fig. 5. The Feynman diagrams of the effective interaction up to third order. The solid and dashed 
lines correspond to the bare Green's function and the interaction, respectively. The diagrams 
enclosed by the dashed line are vertex corrections. The other diagrams are included in RPA. 
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Fig. 8. The calculated as U is varied, at n = 0.53 and for various values of as shown in the 
figure. 
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Fig. 9. The calculated as U is varied, at n = 0.51 and for various values of as shown in the 
figure. 
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Fig. 10. The calculated T^. TOPT and TOPT(RPA-like only) in this figure mean that full diagrams 
and only RPA-like diagrams of anomalous self-energies up to third order are included, respectively. 
See also Fig. 4. 
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Fig. 12. The momentum dependence of the static bare susceptibihty for various tm. The notations 
under the horizontal axis are defined in Fig. 2. These results are obtained for n = 0.55, T = 0.03. 
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Fig. 13. The momentum dependence of the static bare susceptibility for various n. The notations 
under the horizontal axis are defined in Fig. 2. These results are obtained for tm = 0.75, T = 0.02. 
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Fig. 14. The density of states as tm is varied, at n = 0.572, U = 6.0 and T = 0.03. The inset shows 
the details near the Fermi level. 
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Fig. 15. The density of states as n is varied, at tm = 0.75, U = 6.0 and T = 0.02. The inset shows the 
details near the Fermi level. 
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Fig. 16. The real part of the normal self-energy at the Fermi momentum, at tm = 0.75, n = 0.572, 
T = 0.01825 and for various values of U as shown in the figure. 
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Fig. 17. The imaginary part of the normal self-energy at the Fermi momentum, at tm = 0.75, n 
0.572, T = 0.01825 and for various values of U as shown in the figure. 
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